Description of Heavy Quark Systems by means of Energy Dependent Potentials 
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We apply, for the first time, an energy dependent Schrodinger equation to describe static properties 
of heavy quark systems, i.e. charmonium and bottonium. We show that a good description of the 
eigenstates and reasonable values for the widths can be obtained. Values of the radii and of the 
density at the origin are also given. We compare the results to those deduced with a Schrodinger 
equation implemented with potentials used so far. We note that the energy dependence of the 
confining potential provides a natural mechanism for the saturation of the spectra. Our results 
introduce a new class of potentials for the description of heavy quark systems. 
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Wave equations with energy dependent potentials are 
already familiar in physics. The Pauli-Schrodinger equa- 
tion describing a particle in an external electromagnetic 
field [l| or the Hamiltonian formulation of relativistic 
quantum mechanics, in connection with manifestly co- 
variant formalism with constraints U3- Id- LJl - offer two ex- 
amples of energy dependent potentials. Such dependence 
arises from the relativistic nature of the problem. More- 
over, it has been shown that in the context of relativistic 
quantum mechanics, the potentials can be derived from 
fundamental theories, i.e a link exists with field theory 
H . The present work has a different aim, namely to show 
the characteristic features which can be expected from 
energy dependent potentials, when applied to a specific 
physical system. For this purpose, simple potentials ad- 
mitting analytic or semi-analytic solutions are of great 
help since they facilitate the comparison with usual po- 
tentials. 

The introduction of an energy dependence has several 
implications with respect to usual quantum mechanics. 
For example, the conservation of the norm asks for the 
modification of the scalar product |fj. The energy de- 
pendent potentials must fulfill certain conditions in or- 
der to result in a meaningful quantum theory [?]. For 
those potentials which are acceptable, the corresponding 
Schrodinger equation is equivalent to a Schrodinger equa- 
tion including a non-local potential. The non-locality, 
however, is treated more efficiently under the form of an 
energy dependent potential. A detailed analysis of these 
formal aspects is performed in [7] showing that the prop- 
erties of a good quantum theory are well preserved. 

The scope of this letter is to apply, for the first time, 
such an approach to cases of physical interests. A natu- 
ral application is offered by heavy quark systems. Their 



erties of charmonium and bottonium. In particular, we 
present the eigenstates of these systems, the root mean 
square radii, the leptonic widths and the density of the 
S-wave functions at the origin, the latter being a use- 
ful input for quarkonium production in high-energy col- 
lisions. We compare our results to those obtained with 
the potentials known so far to give reasonable descrip- 
tion of charmonium and bottonium properties. We take 
the Buchmiiller and Tye (BT) 8] and the Cornell [3 
potentials as examples. 

For simplicity we consider spherical symmetric poten- 
tials, assuming a linear energy dependence. The corre- 
sponding Schrodinger equation reads (ignoring spin de- 
grees of freedom): 



[-_ V 2 + V{r, E)}^ n4 {r) = E n ^ n4 {r), (1) 



2m 



where the potential used is: 



V(r, £) = (! + jE)f(r) + D P + V 



(2) 



with n, £ denoting the principal and the orbital quan- 
tum numbers respectively and m/2 the reduced mass (we 
work in units fi,=197.3 MeV fm _1 ). To illustrate our aim 
we choose two radial dependences for which analytical 
solutions exist: i) the linear one, i.e. f(r) = Xmr and n) 
the harmonic one, i.e. f(r) — mui 2 r 2 /2. In each case the 
total number of parameters is 5. Note that only four of 
the parameters are free since M(qq) = 2m + Vb- They 
are determined by constraining the energy differences be- 
tween the ground state and the IP, 2S, 2P and 3S states 
to the experimental energies for both systems. The ID 
state is also considered for the cc system. Note that for 
the two chosen radial dependence, only negative values 
of 7 are acceptable Q. 

The term DqP is added to break the degeneracy be- 
tween the states of different angular momentum. Its con- 

tribution is necessary to achieve accurate fits, especially 

in the case of the harmonic oscillator. 
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description has been a quite successful 
potential models in the past 
we show that energy dependent potentials represent a 
new class of potentials that can account for the prop- 
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tained by solving the following quadratic equation: 



El 



2Dl(l + 1) + ahW-f 



(3) 



+ D 2 e 2 (£+ l) 2 -ah 2 uj 2 = 



where D = ^D Q and a = (An + 21 + 3) 2 / 4; while for the 
linear dependence the S-state eigenvalues are obtained 
by solving 

(E n ,o + V ) 3/2 - A n Vn~i\{l + 7 £„,o) = , (4) 

with A n = (-« n ) 3 / 2 -^, a n being the n-th zero of the 
Airy function. 

Eqs. (3) and (4) are non-linear. Therefore one has 
to specify the choice of the eigenvalue in the case when 
several solutions exist. For the harmonic oscillator, only 
positive energies lead to square integrablc wave functions 
(see 0)- Consequently negative solutions are rejected. 
For the linear potential, Eq. (4) can be transformed into 
a cubic equation, which has a single real root in the range 
of parameters used in the present work. 

In the framework of potential models, the leptonic 
widths of the S-states. without radiative and relativis- 
tic corrections |14L h5l Il6| , are given by the van Royen- 
Weisskopf formula |l7j : 



16ire 2 oe 2 



(0)| 5 



(5) 



where R n ,o indicates the radial component of the wave 
function evaluated at the origin, e q the quark charge and 
M(qq) is the physical mass. We have checked that the 
same formula holds when using energy dependent poten- 
tials. 

Table I shows the sensitivity of the results to the varia- 
tion of the parameters, in the case of the harmonic oscil- 
lator radial dependence as an example. The remarkable 
fact is the change in the spectrum pattern brought by 
the energy dependence, which deepens the ground state 
with respect to the other states. The usual harmonic os- 
cillator would be unable to produce the relative positions 
of the IP, 2S, ID and 3S levels, even in the presence of 
the ? term in the Hamiltonian. On the other hand, the 
energy dependent potential with the harmonic oscillator 
radial shape is not able to account for both the ID and 
3S states simultaneously. In this respect, the linear po- 
tential comes out to give much better description of the 
cc and bb properties. 

In Tables II and III we give the results, obtained 
with the linear dependence, for the energies, the leptonic 
widths, the root mean square radius, and the density of 
the IS states at the origin. We compare our results with 
those obtained with previous potentials, in particular the 
BT and the Cornell ones. Their values of the IS energy 
density at the origin are taken from ^^l- We see that a 
very good agreement on the spectrum is obtained. On 
the other hand the width of the IS state as well as the 



SPECTRUM 
IP 

2S 
ID 
2P 
3S 

W.F. PROPERTIES 

Tis 



ris/r2s 

'IS 



PARAMETERS 

lo (fai -1 ) 
Vq (MeV) 
D (MeV) 
■ 1CT 4 (MeV -3 
m(GeV) 



397 397 397 397 

589 589 589 589 

758.5 759.2 706.2 706.4 

839.0 839.6 789.5 789.6 

946.2 947.1 876.7 876.8 



2.6 6.2 2.9 7.7 

1.7 1.5 2.3 1.7 
0.490 0.374 0.469 0.346 



2.7 2.3 3.5 3.0 
0.0 600. 0.0 600. 

2.8 4.3 2.4 3.6 
-4.29 -5.72 -4.75 -6.64 
1.212 1.812 1.138 1.738 



EXP. 

397 
589 
710 
840 
940 



5.26(.37) 
2.4 



TABLE I: Sensitivity of the cc properties to variations 
OF the parameters: The results display the energy differ- 
ences (MeV) with respect to the ground state, the leptonic 
widths (keV), the root mean square radius (fm) in the case 
of the harmonic oscillator radial dependence. Experimental 
data are given for comparison, when available. 



IP 
2S 
ID 
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3S 
Lis 
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This work BT Cornell 



397 
589 
710 
840 
943 
7.51 
2.18 
0.33 
0.761 



424 
601 
716 
872 
1023 
8.0 
2.17 
0.42 
0.810 



428 

590 

713 

871 

1015 

14.33 

2.22 

0.47 

1.454 



397 
589 
710 
840 
940 
5.26(.37) 
2.4 



TABLE II: cc properties: Results obtained with an en- 
ergy dependent potential and the linear radial dependence, in 
comparison with those of Buchmuller and Tye (BT) and 
of the Cornell [jj potentials. The values listed correspond 
to the energies (MeV) with respect to the ground state, to 
the leptonic widths (keV), the root mean square radius (fm), 
the density at the origin (GeV 3 ). The parameters used are 
A=2 foa -1 ,Vb=1940 MeV, D =2 MeV, 7 =-5.007 10" 4 MeV -1 
and m= 1913 MeV. Experimental data are also given. 



ratio of the 2S to the IS widths present the same quanti- 
tative agreement with experiment as previous potentials. 

In this respect, one should remember that the van 
Royen-Weisskopf formula should be multiplied by a fac- 
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BT 


Cornell 


EXP. 


IP 


426 


430 


498 


428 


2S 


563 


560 


591 


563 


ID 


798 


680 


747 




2P 


818 


790 


852 


792 


3S 


895 


890 


936 


895 


4S 


1116 


1160 


1213 


1120 


Tis 


0.93 


1.68 


3.66 


1.32(.05) 


ris/r2s 


2.73 


2.2 


2.17 


2.54 




0.23 


0.23 


0.20 




\Ris(o)\ 2 


3.65 


6.477 


14.08 





TABLE III: bb properties: Same as Table II, but here the re- 
sults correspond to A=0.6 fm"\Vb=2725 MeV, Do = 15 MeV, 
7=-1.06 10" 3 MeV" 1 and m= 5475 MeV. 

tor (l-16a s /37r) corresponding to the (first-order) radia- 
tive corrections 0, • ^ n 01 ^ ^ s shown that af- 
ter the inclusion of the relativistic corrections, potential 
models having |i?i S (0)| 2 =0.8-l (6-8) GeV 3 for the J ftp 
(T) are able to reproduce the leptonic widths. This point 
should be reconsidered in the framework of energy de- 
pendent potentials. Here we shall merely compare the 
zeroth order values with the ones obtained from other 
potentials. In the case of the cc system, the IS energy 
density at the origin is close to the value given by the 
BT potential. On the other hand, for the bb system, it 
is sensibly lower than the BT value. This situation is 
quite strange, since the value of |i?is(0)| 2 is closely re- 
lated to the positions of the S-states, which are very well 
reproduced by the present model (and is not due to a 
peculiar set of parameters). Indeed the density at the 



origin of the lS-state can hardly be increased by 10-15 % 
by choosing a lower coupling constant A. Consequently, 
such a behaviour may point to the necessity of finding a 
more universal radial shape and/or of a more pertinent 
energy dependence. 

In conclusion, we show that a new class of potentials, 
namely energy dependent potentials, give a good descrip- 
tion of charmonium and bottonium properties. We em- 
phasize that the energy dependence is rather naturally 
introduced, since it arises from the relativistic aspect of 
the problem. As a starting point, we have chosen po- 
tentials admitting analytical solutions, at least for the S 
states. This is a great advantage in view of the high non- 
linearity of the differential equation to be solved. How- 
ever, it could be rewarding to study other radial shapes, 
as well as combinations of usual and energy dependent 
potentials. As far as the energy dependence is concerned, 
the linear one has the advantage of simplicity in many re- 
spects as already emphasized in 0. 

We would like to stress a very important feature of the 
energy dependence for confining potentials, namely the 
saturation of the spectra. For usual potential, E n> i in- 
creases regularly toward oo with n or t. In contrast for 
energy dependent confining potentials, E n ^ tends to a 
limit as n and I increase. Thus, the density of states be- 
comes rapidly very large, the splittings getting infinites- 
imally small, which forbids the observation of individual 
states. We believe that this property could have relevant 
implications in many aspects of the quark-model physics. 

Finally, understanding the connection between the 
phenomcnological potentials used here and more funda- 
mental theories of elementary particles deserves further 
studies. The application to other systems is foreseen. 

We thank Yuri Ivanov for providing us with useful infor- 
mation. 
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